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ABSTRACT 


A  digital  computer  program  was  developed  for  performing  the  Bruceton 
statistical  analysis  on  an  IBM  709  or  IBM  7090.  A  brief  description  of  the 
Bruceton  test  is  given  --  followed  by  an  outline  of  the  method  of  analysis 
which  this  program  uses.  Input  and  output  formats  are  discussed. 


INTRODUCTION 


The  purpose  of  this  report  is  to  describe  a  digital  computer  program 
which  was  developed  for  performing  a  statistical  analysis  of  sensitivity 
data  obtained  from  the  Bruceton  Test. 

Part  I  discusses  the  input  and  output  formats.  Part  II  discusses  the 
numerical  calculations  used.  The  appendix  lists  the  program  as  written 
in  Fortran. 

The  Bruceton  Test  --  developed  by  the  Explosives  Research  Laboratory 
at  Bruceton,  Pennsylvania  --  measures  electrical  sensitivity  of  an 
electro-explosive  devi<"-»  at  a  specific  level  of  frequency.  The  lowest 
power  level  should  fall  above  the  "no-fire"  level  for  the  device  being  eval¬ 
uated  and  the  highest  power  level  should  fall  below  the  "all-fire"  power 
level.  This  test  is  preferable  to  other  similar  tests  in  that  it  concentrates 
the  testing  near  the  mean.  Therefore,  in  most  situations,  a  desired 
degree  of  reliability  can  be  obtained  with  a  smaller  sample. 

With  the  data  from  a  Bruceton  Test  the  mean  (m)  or  the  50%  firing 
level  and  the  standard  deviation '((f)  can  be  calculated.  Further  computa¬ 
tion  yields  the  extreme  probabilities  such  as  the  99.9%  and  the  0. 1%  levels. 
These  are  expressed  with  the  confidence  intervals  included. 

This  statistical  analysis  is  based  on  the  Bruceton  Report  (Reference  1). 
It  is  assumed  that  the  calculations  found  in  Part  II  will  be  used  in  conjunc¬ 
tion  with  the  Bruceton  Report  and  the  Franklin  Institute  Memorandum 
(Reference  2)> 
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DISCUSSION 


Part  I  -  Input  and  Output  Formats 


A,  Input  Format 

The  first  ten  data  cards  will  be  the  same  for  any  group  of  runs. 

These  are  constant  values  (G,  H  and  S)  obtained  from  the  graphs  foun  1 
in  the  Bruceton  Report, 

Following  this  set  of  cards,  any  number  of  runs  may  be  included 
provided  each  run  is  in  correct  sequence. 

The  format  for  the  first  data  card  is: 

Spaces  1  through  20  (inclusive)  are  reserved  for  an  identification. 
Spaces  21-25  are  reserved  for  the  number  of  levels  used. 

Spaces  26-35  are  reserved  for  the  common  logarithm  of  the  lowest 

level. 

Spaces  36-45  are  reserved  for  the  constant  interval  between  the 
logarithms  of  the  levels. 

Spaces  46-50  are  reserved  for  the  probability  level  (in  percent). 
Spaces  51-55  are  reserved  for  the  confidence  level  (in  percent). 
Spaces  56-60  are  reserved  for  the  value  of  the  first  (lowest) 
level  (in  milliwatts) . 

Spaces  61-65  are  reserved  for  the  number  of  "fires"  at  that  level. 
Spaces  66-70  are  reserved  for  the  number  of  "no-fires"  at  that 

level. 

Spaces  71-80  are  blank. 

Each  successive  level,  followed  by  the  "fires"  and  "no-fires"  at  that 
level  (in  that  order),  are  placed  on  the  second  card  with  five  spaces  re¬ 
served  for  each  value.  If  more  data  remains  after  the  second  card  is  used; 
a  third,  fourth,  etc,,  may  be  used  with  the  same  format. 

As  a  convenience  those  numbers  above,  which  are  whole  numbers,  need 
not  be  written  with  a  decimal, 

B.  Output  Format 

The  output  will  list  all  the  input  data  as  well  as  every  intermediary 
value  used  in  the  analysis.  The  values  used  in  the  primary  statistics  and 
in  the  secondary  statistics  are  so  labeled. 
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If  P  and  E  equal  the  input  probability  and  confidence  levels  re spec* 
tively»  then  in  the  final  calculations  the  values  (both  in  logarithmic  units 
and  milliwatts) consist,  of 

P%  probability  with  E%  confidence  and 
( 100 -P)%  probability  with  E%  confidence. 
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Part  II  -  Numerical  Calculations 


A.  Primary  Statistics 

Let  N  equal  the  number  of  levels  used. 

Let  X(i),  where  i=0,  1,  2,  ,  ..N,  equal  the  value  of  each  level. 
Let  nx(i)  equal  the  number  of  "fires"  at  each  level. 

Let  n^(i)  equal  the  number  of  "no-fires"  at  each  level. 

Let  c  equal  the  logarithm  of  the  lowest  level. 

Let  d  equal  the  constant. logarithmic  interval  between  levels. 
Let  P%  equal  the  probability  level. 

Let  E%  equal  the  confidence  level. 

-  Nx=fnx(i> 

Ax  *  £  lnx<i> 

Bx  *?  i2nx(i) 


No'?no<l> 
Ao  =  f  ir*o(i> 

Bo  =fi2no<l> 


N  2 


=  c  +  d(  Ax  -  1/2) 
B„ 


Mo  ~  NoBo~An' 


N, 


mQ  =  c  +  d(Ao  - 
B„ 


If  Mj£<  Q.3,  the  smallest  absolute  value  R  = 


mx  -  log  X(i) 
d 


is  evaluated  and  Graph  I  in  the  Bruceton  Report  is  used  to  find  s  . 
rated  in  the  program  (data  cards  3-10)  are  268  points  taken  from  t 
with  M*  =  .05.  .06,  .....  .30  and  R  =  0,  .5,  .1,  .15, . 5 


1/2) 


Incorpo 
is  graph 
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Similarly  sQ  if  mq  <  0. 3. 

the  computation  is  <TX  =  sxd 


<T  =  s  d 
w  o  a 


However,  if  M  £,0.3  then 


0X  S  1.62d  (Mjj  +  0.029) 


01  =  1. 62d(M  +  0 


/xr 

v  Nx-» 


B.  Secondary  Statistics 


m  =  N0m0  + 

N  +  N 
o  x 


^  =  Antilog  m;  m  equals  the  mean  or  50%  firing  level 
and  <5  equals  the  standard  deviation. 


N  =  N0  +  Nx 


n 


when  N  is  an  even  integer 


-  N  +  1 
2 


when  N  is  an  odd  integer 
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C.  Confidence  Interval  (Y) 


Y 


1/2 

^ n  +  *»2)  ^  G2  +  H2kp2^ 


where  k  ,  k  ,  G  and  H  are  obtained  from  the  Bruceton  Report, 

and  kp  are  found  in  the  table  on  Page  19  and  G  and  H  from  Graphs 
Ill,  IV  and  V,  The  table,  in  its  entirety,  is  incorporated  in  the  program, 
and  17  points  are  taken  from  each  of  the  graphs  and  placed  on  data  cards 
1  and  2,  . (S  »  ,2,  ,3,  ,,,,  1,0, 1,5,  2,0,,,.,  5,0),  Straight  line  inter¬ 
polations  are  made  for  values  between  these  points. 


D,  Final  Calculations 

P%  (E%  confidence)  =  m  +  kp  +  Y 
(100-PJ&  (E%  confidence)  =  m—  kp  +  Y 
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APPENDIX 


APPENDIX  A 


FORTRAN  PROGRAM  FOR  A  BRUCETON  ANALYSIS 


BRUCET0N  STATISTICAL  ANALYSIS 


LCGF(W)  *  8«4343*L0GF(W) 

ANTF(W)  *  EXPF  (W/0.4343) 

DIMENSION  Gt50),H(50)  ,X(50> , IDENT120) , FIX (50) , FAX (50) ,Y(50), AB(30, 
115) 

1  READ  INPUT  TAPE  2 , 2, ( G ( I ) , I *2 , 10) , GC 15 ) , G< 20) , G( 25 ) ,G( 30 ) , G(  35) ,G ( 
140) »G ( 45) »G( 50  ) 

2  F0RMAT( 17F4.3) 

3  READ  INPUT  TAPE  2 • 4, ( H (I ) , I *2 , 10) , H( 15 ) , H( 20) ,H( 25 ) , H( 301 , HI  35)  ,M( 
140) »H(45) • H( 50) 

4  F0RMAT ( 17F3.2 ) 

READ  INPUT  TAPE  2 , 7*  ( (AB(J,l),J*5,30),I*l,U) 

7  FORMAT ( 36F2*0/  36F2. 0/36F2. 0/ 36F2 • 0/36F2. 0/36F2.0/36F2.0/34F2. 0) 

5  READ  INPUT  TAPE  2 .6, U DENT ( J ) , J*l ,20) ,N, C »0,P ,E, C X(I) ,FIX( I ) , FAX(  I 
1) ,I*1,N) 

6  F0RMAT(2OA1,I5,2FIO.O,5F5.O,2X/(14F5.O,2X)) 

PI  *  100. -P 

AFIX  *0.0 
AFAX*  0.0 
BF IX*  0.0 


BFAX* 

0.0 

SF  I X* 

0.0 

SFAX* 

1*0 

W*  I 

0.0 

AFIX 

*AF  I X 

♦.  W»FIX(I  +  1) 

AFAX 

*AFAX 

+  W«FAX(I+1) 

BF  I X 

*BF  I X 

+(W**2)*FIX(I+l) 

BFAX 

*BFAX 

*( W*«2)»FAX( I+l) 

SFIX 

*SFIX 

RFIXtI+l) 

SFAX 
1*1  +  1 

*SF  AX 

+.FAXII  +  1) 

PAGE  1 


IFU-N+l)  10,10,11 

11  R*0 

TFIX  *( (SFIX*BFIX)-AFIX**2)/(SFIX**2) 

TFAX* ( ( SFAX*BF AX) -AF  AX**2 )/(SFAX**2) 

13  UFIX*  C+R*  ( AF I X/SFIX-0. 5) 

UFAX*  C+R*(AFAX/SFAX+0.5) 

I F ( TF I X-0 .3)  118,23,23 
118  C I X*TF I X* 100. 

L I X*CI X 

111  WEE  I X*500. 

1*0 

112  WEEIX*MIN1FCWEEIX,ABSF(UFI X-LCGF ( X ( I +1 ) ) ) ) 

i*i+i 

IF ( I-N+l )  112,112,113 
113!  RI  X*ABSF  ( WEEIX/R) 

K!X*XINTF(RIX*22. )+l 
VFIX*(AB(LIX,KIX)/100. ) *R 
24  IF ( TFAX-0. 3)  117,25,25 
117  CAX*TFAX* 100. 

LAX*CAX 

114  WEEAX*500. 

1*0 

115  WEEAX*MINIFIMEEAX,A8SF(UFAX-LCGF(X( I ♦ 1 ) ) ) ) 
IM  +  1 
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BRUCE T0N  STATISTICAL  ANALYSIS 


IFII-N+1)  115,115,116 
i  16  RAX*ABSF ( WEEAX/R) 

KAX*XINTF(RAX*22. )♦! 

VFAX*(ABILAX,KAX)/100. )*R 
G0  T0  26 

23  VF IX*  l.62*R»( TFI X+0.029) *SQRTF ( SFIX/ ( SF IX-i.  ) ) 

S0  T0  24 

25  VF AX*  l. 62*R* ( TFAX+O. 029) *SQRTF ( SFAX/ ( SFAX-1. ) ) 

26  U  *( (SFAX*UFAX)  +  (SFIX*UFIX) ) / ( SFAX+SF I X) 

V  «  SQRTFHSFAX*(  VFAX*»2)  +  SFIX*(  VFIX**2)  )/( SFAX+SFIX) ) 
ZETA*ANTF(U) 

IF ( 99.-P)  51,50*50 

50  KP-99.9-P 
G0  T0  60 

51  KP*{  100«-P)*100«< 

60  IF ( 99.-E )  53,52*52 

52  KX*99. 9-E 
G0  T0  61 

53  KX*( 100.-E )*100» 

61  Y  ( 49) *0. 0 

V ( 24 ) *0.675 
Y(9)*1.282 
Y ( 4) =1 • 645 
Y(0)*2.326 
Y ( 10 ) *3.090 

Y(l)*3.719  PAGE  2 

S=  V/R  f 

Z  *  S*10. 

SX*  SFIX+SFAX 
SX2  *(SX+i.)/2. 

KS2  *  SX2 
SX3  *  KS2 
IFU-2.)  100*14,14 

14  IFCZ-50. 115,15,100 

15  IFIZ-2. 1100416, 16 

16  IF(Z-I0.)40»i7»18 

17  K*Z 
A*G(K) 

B*H( K) 

G0  T0  19 

40  K*Z 
S2*K 

A* t G t K+ l )-G( K)  )*(Z-S2)  ♦  GIK) 

B*(H(K+l)-H(K)  )*IZ-S2)  ♦  HIK) 

G0  T0  19 

18  SI*  Z/5. 

XI*  SI 
K*5*Kl 
S2*K 

A  *((GIK*5)-GCK))/5.)*(Z-S2)^  GIX) 

B  *< (H(K+5)-H(K))/5.)*(Z-S2)+  HIK) 

19  T  *YCKX)*USX3+l.2)/SX3)*SQRTF(  I  A*»2«- (  B*»2)  MYCKP  )**2) ) /SX3)*V 
CFll  *U+Y ( KP) * V+T 

CFL2  *U-Y ( KP )  • V~T 

CFN1  *  ANTF(CFLl)  .9. 

CFN2  *ANTF (CFL2) 


BRUCET0N  STATISTICAL  ANALYSIS 


WRITE  OUTPUT  TAPE  3, 20 , (I  DENT ( J ) , J*  1 , 20) , P , E, ( XII) ,F I X( 1 1 , FAX (I) , I 
1*1, N) 

20  F0RMATIIHI,IOX,2OA2,5X,13HPR0B.  LEV.  *  F7.3,5X,13HC0NF.  LEVi  *  F7. 
13//10X, ilHTEST  LEVELS, 5X,12HN0. OF  F I  RE S,  5X , 1 1HN0 . 0F  FAILS//! 10X,F7 
2.3, IOX,F7.3,10X,F7.3) ) 

WRITE  OUTPUT  TAPE  3 , 21 , SF AX, SF I  X, C, R, AFAX , AF I X,BFAX, BF I Xi TFAX , TF I X 

I, UFAX,UFIX,VFAX,VFIX,U,V,ZETA,P,Pl,E,S,SX,SX3,Y!KP),Y(KX)  ,A,B,T,.P, 
2E,CFLI,CFNi»Pl ,E,CFL2,CFN2 

21  FORMAT! 1H0, 10X ,9H  N0,NX  *  2F11.5//26H  SPECIAL  PARAMETERS  C,D  *  2FI 

II. 5/66H  PRIMARY  STATISTICS  AO , AX, B0, BX, MO , MX, M01 , MX1 , SIGMAO, SI 

2GMAX  -  ///5!13X,Fll45,5X,Fll.5/),10X///45H  SEC0NOARY  STATISTICS  ME 
3AN, STAND. DEV., ZETA  *  3F11.5//16H  PR0B.  LEVELS  *  2FI1.5, 14HC0NF.  LE 
4VEL  *  FI 1 • 5/5X  »I9HS»N» N1 » KP» KX*G» H  *  /I0X,6(F  11.5,3X)/10X,F11.5,3X 
5/25H  CONFIDENCE  INTERVAL  Y  *  F1I.5//1 IX,F7.3,5H  WITH,F7.3,8H  CONF 
6*  FII.5,9HL0G  UNITS,4H  *  F11.5/)) 

GO  TO  102 

100  WRITE  OUTPUT  TAPE  3, 90, ! I  DENT ( J ) , J*I • 20) , P , E, I  X! I) ,F I X! I ) ,FAX ( I ) , I 

1*1, N) 

90  FORMAT! 1H1,10X,20A2,5X, 13HPR0B.  LEV.  *  F7.3t5X,13HC0NF.  LEV*  *  F7. 
13//10X, IlHTEST  LEVELS, 5X, 12HN0. OF  FIRES, 5X, 1 IHNO. OF  FAILS//! 10X,F7 
2.3, 10X,F7.3*I0X,F7.3) ) 

WRITE  OUTPUT  TAPE  3,91,S 

91  FORMAT! 1H0, 10X ,44HTHI S  ANALYSIS  WAS  NOT  COMPLETED  BECAUSE  S  *  Fll. 
15,i03H  IT  IS  OUT  OF  THE  RANGE  FROM  .2  TO  5  AND  HENCE  G  AND  H  CANNO 
2T  BE  DETERMINED^  CHECK  LEVEL  DISTRIBUTION. I 

102  SO  TO  5 

END! 1,1, 0,0, 0,0, l, l, 0,0, 0,0, 0,0,0) 
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